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Abstract 

Background: In ovo electroporation is a widely used technique to study gene function in developmental biology. 
Despite the widespread acceptance of this technique, no genome-wide analysis of the effects of in ovo 
electroporation, principally the current applied across the tissue and exogenous vector DNA introduced, on 
endogenous gene expression has been undertaken. Here, the effects of electric current and expression of a GFP- 
containing construct, via electroporation into the midbrain of Hamburger-Hamilton stage 10 chicken embryos, are 
analysed by microarray. 

Results: Both current alone and in combination with exogenous DNA expression have a small but reproducible 
effect on endogenous gene expression, changing the expression of the genes represented on the array by less 
than 0.1% (current) and less than 0.5% (current + DNA), respectively. The subset of genes regulated by electric 
current and exogenous DNA span a disparate set of cellular functions. However, no genes involved in the regional 
identity were affected. In sharp contrast to this, electroporation of a known transcription factor, DmrtS, caused a 
much greater change in gene expression. 

Conclusions: These findings represent the first systematic genome-wide analysis of the effects of in ovo 
electroporation on gene expression during embryonic development. The analysis reveals that this process has 
minimal impact on the genetic basis of cell fate specification. Thus, the study demonstrates the validity of the in 
ovo electroporation technique to study gene function and expression during development. Furthermore, the data 
presented here can be used as a resource to refine the set of transcriptional responders in future in ovo 
electroporation studies of specific gene function. 



Background 

Determining the function of genes involved in embryo- 
nic development requires manipulation of gene expres- 
sion in a spatially and temporally restricted manner. 
Whilst transgenesis in mice is used to this end, it is a 
time consuming and costly technique. Moreover, achiev- 
ing the required spatiotemporal targeting via conditional 
transgenesis is not always achievable [1]. In contrast, the 
highly precise localisation and timing of expression con- 
struct insertion made possible by the in ovo electropora- 
tion technique provides an economical and efficient 
alternative to transgenesis [2]. These advantages make 
in ovo electroporation a widely used technique to study 
gene function during development. 
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In ovo electroporation is employed in both gain and 
loss of function studies and reporter studies to analyse 
gene function and regulation during development. This 
method has been applied to all areas of developmental 
biology, including neurogenesis and neural differentia- 
tion [3,4], axon outgrowth and guidance [5] axonal pat- 
terning in the developing Umb [6], somitogenesis [7,8], 
skeletal muscle development [9] and eye development 
[10,11]. 

In ovo electroporation has been used to study develop- 
ment in mouse embryos [12] along with other verte- 
brates, including Danio rerio [13] and Xenopus laevis 
[13,14] and non-vertebrates such as Drosophila melano- 
gaster [15] and Ascidiacea [16] (for reviews, see [17,18]). 
However, in ovo electroporation has been most widely 
applied to avian embryos because the avian embryo 
develops in ovo and has a planar topology. These char- 
acteristics greatly facilitate injection of the DNA 
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construct, placement of electrodes and incubation of the 
electroporated embryo. 

Independent of the experimental organism, electro- 
poration involves microinjecting a gene expression 
construct, typically into a natural body cavity such as 
the neural tube. Subsequently, electrodes are placed 
flanking the site of injection and an electric current is 
applied across the embryo in the form of a rapid series 
of square wave pluses. This electric field transiently 
disrupts the stability of the plasma membrane, creating 
pores in the cell membrane. The negatively charged 
DNA constructs migrate towards the positive electrode 
and enter the cells in their path via these pores. The 
tissue adjacent to the negatively charged electrode 
remains untransfected, providing an internal control 
[19]. This method enables defined tissues to be tar- 
geted, by location of DNA injection and positioning of 
electrodes, at precise times of development. For more 
detailed information about specific parameters required 
refer to [2,20]. Further advances now enable focal elec- 
troporation and electroporation of different constructs 
in close proximity, increasing the precision and com- 
plexity of in ovo electroporation studies. This techni- 
que uses beads soaked in the DNA construct, which 
are microsurgically implanted into the embryo instead 
of microinjection [21]. 

In ovo electroporation is most commonly applied to 
Hamburger-Hamilton (HH) stage (st) 10 to 20 embryos. 
Embryos older than HH st20 have more compact tis- 
sues and increased tissue layers, making microinjection 
and in ovo electroporation more difficult. For older 
embryos ex ovo explant electroporation is an alternative 
method [22]. 

Typically, the construct used for in ovo electropora- 
tion contains the gene of choice along with green fluor- 
escent protein (GFP), or another marker, in order to 
identify the cells that have taken up the construct. A 
popular construct to obtain transient expression is 
pCAP'IRES-GFP containing the P-actin promoter, a 
cytomegalovirus (CMV) enhancer, a polylinker for 
inserting the desired gene followed by an internal ribo- 
somal entry site {IRES) and GFP [23]. Expression of 
the translated product of the construct - for example, 
GFP - can be detected 2.5 hours after electroporation 
and peaks around 20 to 24 hours [24]. Expression of 
transient constructs can be maintained for 3 to 11 
days [2,25]. Constitutive expression can be obtained by 
integrating plasmids into the genome using methods 
such as transposon-mediated gene transfer [26]. Using 
constructs with inducible promoters - for example, the 
tetracycline on tetracycline off system - enables further 
control over the timing of exogenous DNA expression 
[27]. It is also possible to use cell-type-specific enhan- 
cers [2]. 



In addition to overexpression studies, in ovo electro- 
poration can also be used to carry out loss-of-function 
studies using RNA interference [28,29] or dominant 
negative constructs [30]. Constructs containing the gene 
of interest linked to either a repressor or an activator 
enable investigation of the transcriptional activity of 
genes in vivo [31]. As well as investigating gene func- 
tion, in ovo electroporation can be used to study the 
activity of promoters during development in vivo using 
reporter constructs [32,33]. 

Traditionally, the downstream effects of in ovo electro- 
poration of a gene of interest have been analysed using 
in situ hybridisation and antibody staining of a few 
select gene products. Parallel analysis of these genes in 
embryos electroporated with GFP alone is used as a 
control to identify any non-specific effects. This 
approach has the major disadvantage of surveying only a 
limited number of genes for qualitative transcriptional 
changes. To date, no effect of the in ovo electroporation 
technique on gene expression has been reported. How- 
ever, evidence from other techniques, such as electro- 
chemotherapy and electrogene therapy on malignant 
melanoma cells, indicates that both electric current and 
exogenous expression of DNA have a small but signifi- 
cant affect on endogenous gene expression [34]. 

The other potentially influential component of in ovo 
electroporation is the introduction of DNA constructs, 
which results in significant levels of exogenous DNA 
within the electroporated cells. The cellular response to 
expression of plasmid DNA has been analysed in several 
cell lines: Chinese hamster ovary epithelial cells (CHO- 
Kl), mouse fibroblast cells (NIH3T3), human embryonic 
kidney cells (HEK293) and several melanoma lines [35]. 
Exogenous DNA expression in these cell lines was 
shown to have differing affects, triggering a DNA 
damage response in CHO-Kl and NIH3T3 cells but 
causing no cellular response in HEK293 cells or mela- 
noma cell lines [35]. These differences in response to 
exogenous DNA expression are postulated to be due to 
species differences [35]. How chicken cells respond to 
exogenous DNA expression or indeed how multicellular 
living organisms respond is yet to be examined. 

Despite these reports suggesting exposure to current 
and exogenous DNA within cells can cause changes in 
endogenous gene expression, there has been no sys- 
tematic analysis of the effects of either of these compo- 
nents within the developing embryo. Understanding the 
experimentally induced changes caused by in ovo elec- 
troporation has recently become extremely pertinent 
since the sequencing of the chicken genome and the 
availability of a highly representative chicken genome 
microarray (Affymetrix GeneChip Chicken Genome 
Array) now enables the coupling of in ovo electropora- 
tion to genome-wide analysis. This strategy provides an 
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opportunity to use in ovo electroporation to unravel the 
role of genes during the embryonic development in a 
high-throughput and genome-wide manner. If this strat- 
egy is to be used, a better understanding of the effects 
of this technique are required. To address this, we car- 
ried out microarray analysis of the effects of both the 
electric current and exogenous DNA expression asso- 
ciated with in ovo electroporation. 

Results 

Experimental strategy 

To investigate the effects of the electric current and 
exogenous DNA expression associated with in ovo elec- 
troporation on endogenous gene expression of the 
developing chicken embryo, we used microarray analysis 
to compare the gene expression profile of four pooled 
samples of ventral lateral midbrain (VLM) tissue. The 
four samples were unelectroporated VLM (VLM), VLM 
exposed to electric current only (VLMi), VLM exposed 
to expression of GFP and electric current (VLMg) and 
VLM exposed to exogenous expression of a regulatory 
gene, DmrtS {Doublesex and Mab related transcription 
factor like 5) and electric current (VLMd) (Figure 1). 
The current was generated with an ECM830 Electro 
Square Porator (5 x 50-ms pulses of 12 V). The GFP 
construct used was pCA^-IRES-GFP, Electroporation 
was carried out at HH stlO, and embryos were collected 
24 hours later. Thus, our experiment was specifically 
designed to interrogate the set of genes 'stably regulated 
following electroporation, rather than those responding 
to the initial treatment. We reasoned that genes identi- 
fied at this 24 hour time point would be more likely to 
have a significant contribution to altering long-term cell 
identity as opposed to those genes whose expression 
changed rapidly in the minutes following electrical 
insult. 

Effect of in ovo electroporation on endogenous gene 
expression 

The distribution of gene expression values across the 
whole microarray was monitored before and after nor- 
malisation by quantile grouping (box plots) and was 
found to be similar for each set of biological replicates 
(data not shown). The variation in biological replicates 
was also analysed by principal component analysis 
(PCA), which measures the variation in expression levels 
between the microarrays. PCA showed all biological 
replicates cluster together, demonstrating that the biolo- 
gical replicate datasets are similar and reproducible. As 
well as biological replicates clustering together, three 
experimental conditions, VLM, VLMi and VLMg, cluster 
together, indicating there is little difference in gene 
expression between these three conditions. In contrast, 
exogenous expression of the regulatory gene DmrtS 



caused the VLMd samples to cluster separately from the 
other conditions (Figure 2). This preliminary analysis of 
global gene expression levels indicates that exposure to 
current or current + GFP has little effect on endogenous 
gene expression, whilst addition of a regulatory gene 
causes a much larger change in endogenous gene 
expression. 

Analysis of the array data to identify genes that show 
significant differential expression {P > 0.05, >2. 0-fold 
change) upon exposure to current, current + GFP or 
current + DmrtS, when compared to the unelectropo- 
rated tissue, identified that only 41 of the 32,773 oligo- 
nucleotides spotted on the array were changed by 
exposure to current and 176 by current + GFP, This 
equates to a change in gene expression of 0.1% (current) 
and 0.5% (current + GFP), of the oligonucelotides found 
on the array, respectively. By comparison, exogenous 
expression of a regulatory gene, DmrtS, led to a signifi- 
cant change {P > 0.05, >2.0-fold change) in the expres- 
sion of 479 oligonucleotides or 1.5% of the 
oligonucleotides on the array. For further analysis, these 
oligonucleotide lists were processed to remove probe set 
duplicates and oligonucelotides that did not correspond 
to annotated genes (Table 1). 

Effect of current on endogenous gene expression 

Only 41 of the 32,773 oligonucelotides were significantly 
differentially expressed 24 hours after exposure to cur- 
rent alone. This corresponds to 21 annotated endogen- 
ous genes (Figure 3). There was an equal proportion of 
up- and down-regulated genes (10 and 11, respectively) 
and the range of their fold changes was -3.8 to 3.0 
(Table 2). Thus, the application of an electric current 
alone to developing neuronal cells within the ventral 
midbrain has a small effect on the constituents of the 
cell transcriptome. These findings are consistent with 
other in vitro studies that suggest current alone has a 
minimal impact on cellular identity and specification 
[34]. 

We next looked at the composition of the 21 respon- 
ders in this group to determine the set of cellular activ- 
ities affected. The 21 responders belong to a disparate 
array of biological functions (Table 3). The small num- 
ber of genes meant it was not possible to identify biolo- 
gical processes that were significantly enriched in this 
population using Gene Ontology (GO) analysis. 

Upregulation of Heat shock protein 2S {Hsp2S) was 
observed. Two members of the ubiquitination pathway, 
uber2r and ubap2, also responded upon exposure to 
current. Five of the 21 responders are involved in 
response to oxidative stress. However, no major oxida- 
tive stress components show a change in expression. 
Four of the 21 genes are reported to respond to pH/ 
redox changes. Five responders are involved in 
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B. Tissue processing and microarray 



RNA extracted from VLM 
and converted to cDNA 



Samples hybridised 
to array 



WT 



Current 



Current + GFP 



Current + Dmrt5 



Data processing and 
Statistical analysis 



1 . Mas 5 

2. Normalisation _ 

3. Filter 

4. Statistical analysis 



List of genes differentially 
expressed between: 



1. WT and Current 
^2.WT and Current + GFP 
3. WT and Current + Dmrtb 



Figure 1 Experimental strategy. (A) Dissection of the ventral lateral midbrain (VLM) region, (i) In situ sagital view of a HH stl6 electroporated 
embryo expressing pCAp-lRES-DmrtS construct (dorsal view in the inset). Red lines indicate the midbrain (Mb) region, which was dissected out. 
(ii-iv) Coronal sections of midbrain. White lines mark the VLM region. This region was isolated from (ii) control embryos (VLM), (iii) VLM exposed 
to current (VLMi), (iv) VLM exposed to current + GFP (VLMg), and VLM exposed to current + Dmn5 (VLMd; image not shown) for investigation of 
transcriptional profiles by microarray analysis. (B) Tissue processing and microarray analysis. Six VLM tissues were pooled for each biological 
replicate, and three biological replicates were used for each condition. cDNA was isolated from these pools and hybridised to the Affymetrix 
Chicken Genome Array. Following MASS, normalisation and filtering, genes whose expression differed significantly between the wild type (W; 
VLM) and VLMi, VLMg and VLMd were identified by one-way ANOVA. 



regulation of apoptosis, either to promote (two genes) 
or inhibit (three genes) this process. For all three of 
these biological functions, regulators of processes 
rather than determiners are affected, suggesting that 
the effect of the electric current on biological pro- 
cesses is minimal. 



Effect of current and vector DNA on endogenous gene 
expression 

The set of genes derived from the control tissues (VLM) 
was compared with their electroporated counterparts 
expressing GFP (VLMg). This analysis revealed that only 
176 of the 32,773 oligonucleotides represented were 
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X axis = Principal component 1 (greatest variation) 
Y axis = Principal component 2 
Z axis = Principal component 3 

Figure 2 Principal component analysis to identify genome- 
wide transcriptional variation caused by in ovo electroporation. 

This plot shows the variation between the samples; each dot 
represents the global gene expression of a single microarray. The 
greatest variation is measured on the x-axis, then the y- and z-axes, 
respectively. Wild-type VLM, VLMi and VLMg samples cluster 
together on the x-axis, indicating that there is little variation in the 
gene expression between these three conditions. VLMd, in which 
the regulatory gene DmrtS is exogenously expressed, clusters 
separately, indicating larger genome-wide variation between this 
sample and the others. 



significantly differentially expressed {P > 0.05, >2.0-fold 
change) MAhich corresponds to less than 0.5% of the 
transcriptome. This corresponds to 111 endogenous 
genes (Table 4). Of these 111 responders to current + 
GFP, only 10 were upregulated, Mrith the remaining 101 
being downregulated. All of the upregulated genes 
showed a 2.0- to 3.3-fold change except for BSFB 
{beaded filament structural proteinl, filensin), which 
was upregulated 7.8-fold. BSFB is a filament protein and 
forms a structural constituent of eye lens [36], and the 
significance of its upregulation is not known. All down- 
regulated genes showed between a 2.0- and 3.6-fold 
change in expression. Thus, the application of an elec- 
tric current in combination with expression of 



exogenous DNA in the form of GFP has a greater effect 
on neuronal cells than current alone. However, this is 
still a much smaller effect than exogenous expression of 
a regulatory gene, such as DmrtS, which causes a change 
in gene expression of 309 endogenous genes. 

We next looked at the composition of the 111 respon- 
ders in this group to determine the cellular activities 
affected. GO Tree Machine (GOTM) analysis was used 
to identify biological processes that were significantly 
enriched in this population. The biological function 'cel- 
lular component organisation' was significantly affected 
by exposure to current + GFP, Within this category reg- 
ulation of cellular component organisation', 'organelle 
organisation' and 'cytoskeletal organisation' were signifi- 
cantly affected. These changes in organelle organisation' 
included the subterms 'chromatin organisation' and 
'chromatin modification', which were significantly over- 
represented in the tissue exposed to current + GFP, 
This indicates that current + GFP may affect the cellular 
organisation (Figure 4A). This finding was also identified 
by Ingenuity Pathway Analysis (IPA). Using IPA, the 
most significant biological functions within the list of 
responders are 'cell-cell signalling and interaction', 'cel- 
lular assembly and organisation' and 'cellular function 
and maintenance' (Figure 4B). 'Cell death' is the ninth 
most significant pathway, with 9 out of the 111 genes 
found to be involved in this process. However, the 
changes in these apoptotic genes include ones that 
could both promote and reduce apoptosis (Figure 4B). 

To identify if any of the changes in gene expression 
upon exposure to current + GFP elicit a toxic affect, 
such as oxidative stress or apoptosis, IPA toxicity analy- 
sis was used (Figure 5). This analysis found some genes 
involved in pro-apoptosis and p53 signalling in the set 
of 111 responders. However, the number found is less 
than would be expected to occur in a random set of 111 
genes, suggesting this is not statistically significant. 
Thus, in ovo electroporation does not induce toxic 
effects within the embryo as defined by the pathway 
analysis software. 

In ovo electroporation does not affect the regional 
identity of the VLM 

Given the high number of patterning marker genes that 
are well described for the VLM, we were able to analyse 



Table 1 Differential expression of genes upon exposure to current, current + GFP and current + DmrtS 



Condition Number of oligonucleotides Percentage of oligonucleotides showing Number of annotated endogenous genes 
differentially expressed differential expression differentially expressed 



Current 


41 


0.13 


21 


Current + 


176 


0.54 


111 


GFP 








Current + 


479 


1.46 


309 


DmrtS 
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Figure 3 Genes showing differential expression in response to 
current. Genes showing differential expression wlien exposed to 
current and tineir fold change. 



Table 3 Biological functions of genes that responded to 
current 



the effects of current and expression of GFP on pattern- 
ing of this region. The VLM is lateral to the floor plate, 
in the Nkx6.2 expressing region; therefore, we expected 
markers such as Nkx6,2 to be present in the microarray. 
Analysis of the transcriptional profile obtained for the 
VLM was reflective of the known genetic profile of the 
region and the expression levels of the ten marker genes 
were similar in the VLM, VLMi and VLMg samples 
(Figure 6A). 

Indeed, statistical analysis using one-way analysis of 
variance (ANOVA) shows that these marker genes are 
not differentially expressed between the three samples 
to a P-value of 0.05. These data are illustrated in Figure 
6B, which is a box plot of all ten marker genes. All the 



Table 2 Genes showing differential expression upon 
exposure to current and their fold change 

Gene name 



SPIN! 

HNRNPK 

CHDl 

TXNLl 

UBE2R2 

USP9X 

UBAP2 

SYNl 

MIER3 

R0B02 

SGMSl 

CAW 

TRIM36 

HLF 

ACCNl 

ATPJ5 

FST 

MA21L2 
HSP25 
GPR137B 
TMEM118 



Fold change 

-3.8 
-3.0 
-3.4 
-3.2 
-2.9 
-2.6 
-2.4 
-2.6 
-2.4 
-2.3 
-2.3 
2.3 
2.0 
2.0 
2.0 
2.1 
2.3 
2.4 
2.5 
2.6 
3.0 



Biological category 



Genes 



Response to oxidative stress 
Response to pH/redox changes 



Hnrnpk, Usp9x, Hit, Txnll, Hsp25 
CoW, Atp5j, Accnl, Txnll 



Regulation of apoptosis 


Accnl, Hlf, Usp9x, Hnrnpk, Sgmsl 


Table 4 Genes showing differential expression upon 
exposure to current + GFP and their fold change 


Gene name 


Fold change 


MYHW 


-2.0 


RCBTBl 


-2.0 


BAZIB 


-2.0 


ZFHX3 


-2.0 


EP400 


-2.0 


GSPTl 


-2.0 


YLPMl 


-2.0 


ANKHD1-EIF4EBP3 


-2.0 


SGMSl 


-2.0 


MLL 


-2.1 


WAC 


-2.1 


PITPNM3 


-2.1 


TNRC6B 


-2.1 


SMPD3 


-2.1 


EPN2 


-2.1 


CCDC88A 


-2.1 


LOC427360 


-2.1 


TLN2 


-2.1 


DGCR8 


-2.1 


L0C4 18437 


-2.1 


TMTC3 


-2.1 


NIPBL 


-2.1 


PTCHDl 


-2.1 


BEND7 


-2.1 


GNPTAB 


-2.1 


SENP5 


-2.1 


NKTR 


-2.1 


FBXW8 


-2.1 


NISCH 


-2.1 


ACPL2 


-2.1 


LOC429468 


-2.2 


KCNK5 


-2.2 


TBLIX 


-2.2 


UBE4B 


-2.2 


NTRK3 


-2.2 


MLL3 


-2.2 


ALS2 


-2.2 


PHIP 


-2.2 


PPP6R3 


-2.2 


SPEN 


-2.2 


ERCl 


-2.2 
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Table 4 Genes showing differential expression upon 
exposure to current + GFP and their fold change 



(Continued) 


ABCB 1 1 


-2.2 


PRPF38B 


-2.2 


LRRNl 


-2.2 


IL17RD 


-2.2 


TXNLl 


-2.3 


CHDl 


-2.2 


USP40 


-2.3 


PI XNR? 


-2.3 


Ific rRNA 

1 UJ 1 1 \ l V/\ 


-2.3 


SPTBNl 


-2.3 


PCDH19 


-2.3 


PDCD6 


-2.3 


FZD7 


-2.3 


RASGEFIA 


-2.3 


NAV2 


-2.3 


EP300 


-2.3 


GLS 


-2.4 


ATRX 


-2.4 


FAM130A2 


-2.4 


MDNl 


-2.4 


RPRDIA 


-2.4 


VCAN 


-2.4 


ARHGAP18 


-2.4 


JAGl 


-2.4 


BAZ2B 


-2.4 


HNRNPK 

/ // V/ 1/ V/ /\ 


-2.5 


CKAP5 


-2.5 


SYNE2 


-2.5 


ASPM 


-2.5 


PLA2G2E 


-2.6 


FIGN 


-2.6 


HERC2 


-2.6 




-2.6 


5ALL3 


-2.6 


TRAF3 


-2.6 


AKAP13 


-2.6 


RRM2S 


-2.7 


/ OC 7 7 692 7 


-2.7 


SON 


-2.7 


A DUl^ AD 7 V 


-2.7 


AKAP9 


-2.8 


VGLL4 


-2.8 


ZEB2 


-2.8 


LRP2 


-2.8 


Ube2r2 


-2.9 


UBAP2 


-2.9 


YTHDCl 


-3.1 


TRIO 


-3.1 


GUCYIA2 


-3.3 


LM07 


-3.3 



Table 4 Genes showing differential expression upon 
exposure to current + GFP and their fold change 

(Continued) 



PDZRN3 


-3.4 


SFRSI8 


-3.4 


HMGA2 


-3.5 


PRTG 


-3.5 


FH0D3 


-3.5 


VOL 


-3.6 


SPIN 


-3.6 


DLGI 


-3.7 


CACNG3 


2.0 


DNTTIPI 


2.0 


HUSI 


2.0 


PSPH 


2.5 


KCTDC7 


2.8 


WNT3A 


2.8 


EIF4AI 


3.1 


EGFL7 


2.9 


CIRBP 


3.3 


BFSPI 


7.8 



box plots have almost the same mean value and all 
overlap (Figure 6B). These data demonstrate that neither 
electric current nor current -h GFP affect the regional 
identity of the tissue. 

Discussion 

In ovo electroporation is an extremely powerful techni- 
que to investigate the function of genes and regulatory 
regions during development. With advances in methods 
used to analyse electroporated embryos, such as micro- 
array analysis and next generation sequencing, a better 
understanding of the effects of this technique are neces- 
sary. To this end we carried out a genome-wide analysis 
of the effect of this technique on endogenous gene 
expression. 

Our analysis has established that the electric current 
used during electroporation (5 x 50-ms pulses of 12 V) 
has a minimal affect on gene expression, causing a 
change in expression of only 21 endogenous genes. The 
upregulation of Hsp25, which encodes a heat shock pro- 
tein known to be upregulated 24 hours after heat shock 
in response to protein aggregation [37], combined with 
changes in two members of the ubiquitination pathway 
indicate that some protein denaturation may have 
occurred upon exposure to current. 

We have also established that the current and expres- 
sion of exogenous DNA, in the form of GFP, has a 
small but statistically significant effect on gene expres- 
sion, with no toxicity pathways significantly activated. 
The expression of exogenous DNA has a greater affect 
on endogenous gene expression than exposure to 



Farley et al. Neural Development 201 1, 6:17 
http://www.neuraldevelopment.eom/content/6/1 /1 7 



Page 8 of 1 2 



3.0 

.2.5 



^ 2.0 



CO 

> 



1.5 



o 



1.0 
0.5 



0.0 



■ Expected 
□ Observed 




Iling 


on 


mbly 


ation 


notion 


ance 


ment 


CO 




0 


CO 


c 


Q_ 


ign 


rac 


r Ass 


ani 


Fur 


nte 


elo| 




-1— » 




ular 


"co 


Dev 






_co 


O 




ell-C 
and 


CellL 


and 


Cell 


and 


ular 


O 












Cell 




Q. 
X 
LU 

0 
c 

0 

CD 



O 2000-2010 Ingenuity Systems. Inc. All rights reserved. 



0 
O 
>^ 

O 
0 

o 



0 

o 



E 

CO 

o 

CO 
-I— » 

0 

■g 
"o 
< 
o 
c 

"E 
< 




CO 

^ "CO 

= ^ 
0^ CO 

0 

o 



Figure 4 Effect of electric current + GFP expression on endogenous gene expression of the VLM. (A) Genes showing differential 
expression when exposed to current + GFP and their fold changes. For gene names refer to Table 3. (B) GOTM analysis showing biological 
function and expected and observed number of genes in each category compared to the number expected from a random set of 1 1 1 genes. 
(C) IPA functional analysis showing the top ten biological functions enriched in the VLMg compared to VLM samples. 



current alone, with 111 genes affected by exposure to 
current + GFP, 

Interestingly, exogenous DNA expression leads to the 
downregulation of 101 genes, in contrast to the upregu- 
lation of 10 genes, suggesting that exogenous DNA 
expression, in these conditions, results in a repression of 
endogenous gene expression. 

There has been surprisingly little investigation into 
the cellular response to exogenous DNA expression. 



Husl, a gene involved in the ataxia telangiectasia- 
mutated (ATM) and ATM and Rad3-related signalling 
network, has been shown to be upregulated in 
response to exogenous DNA expression in CHO-Kl 
and NI3T3 cells [35] and in response to current + GFP 
in this experiment. IPA and GOTM analysis indicate 
that the exogenous expression of DNA causes changes 
in 'cellular organisation functions, including 'cytoskele- 
tal organisation' and 'chromatin organisation'. Despite 



Farley et al. Neural Development 201 1, 6:17 
http://www.neuraldevelopment.eom/content/6/1 /1 7 



Page 9 of 1 2 



Percentage 



G2/M Transition of the Cell Cycle 
NFkB Signaling Pathway 
Pro-Apoptosis 
PPARc/RXR Activation 
Hypoxia-Inducible Factor Signaling 
VDR/RXR Activation 
TR/RXR Activation 
TGF-p Signaling 
p53 Signaling 



Mechanism of Gene Regulation by 
Peroxisome Proliferators via PPARa 



Aryl Hydrocarbon Receptor Signaling 



RAR Activation 



Oxidative Stress Response Mediated by 
Nrf2 



o o o 



-log(p-value) 



Figure 5 Ingenuity Pathway Analysis toxicity analysis. The 

genes showing differential expression upon exposure to current + 
GFP are not significantly involved in any type of toxicity. The bar 
chart shows the percentage of genes involved in each type of 
toxicity (left axis). Numbers above the bars are the number of genes 
that would equate to 100% for each type of toxicity. The orange 
points show the -log(P-value) (right axis) for each category, and thus 
indicate whether a significant number of genes are found within a 
category to suggest an actual toxicity affect. A -log(P-value) of 1.3 
would be considered significant; all -log(P-values) are lower than 1.3 
and therefore not significant. 



the downregulation of several genes involved in 'chro- 
matin organisation and modification', this does not 
appear to have had an appreciable effect on gene 
expression since only 0.5% of the genes show a change 
in gene expression after exposure to current + GFP. 

The absence of any effect on the regional identity of 
the VLM indicates that any changes seen were not suffi- 
cient to affect the specification of the tissue at the time 
point investigated here. 
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Figure 6 Regional identity of VLM unaffected by current alone 
or current + GFP. (A) Expression profile of ten marker genes of the 
VLM obtained from the microarray analysis of VLM (wild type (V\^), 
VLMi and VLMg. In all conditions, marker genes are seen at the 
expected levels, and there is little difference between the three 
conditions. (B) Box plot showing the normalised expression values 
for all ten marker genes. All conditions have the same median and 
error bars overlap, indicating that the ten marker genes are not 
significantly differentially expressed between the three conditions. 
This shows the regional identity is not significantly affected by 
exposure to current or current + GFP. Statistical analysis using one- 
way ANOVA also shows that these marker genes are not 
differentially expressed. 



In contrast to the above findings, current + DmrtS has 
a far larger affect on gene expression, causing a change 
in expression of 309 genes. These changes in gene 
expression are associated with cell specification, a 
known function of DmrtS (unpublished data). Compari- 
son of the range of fold changes seen in genes affected 
by current + DmrtS (-15 to +1,125), current (-3.8 to 
+3.0) or current + GFP (-3.7 to +7.8) further highlight 
the minimal affect that the current and exogenous DNA 
expression used in this technique have on endogenous 
gene expression when compared to the affect seen when 
a regulatory gene is exogenously expressed. 

Three genes are affected in all experimental condi- 
tions, current, current + GFP and current + DmrtS 
(Figure 7). These are Hnrnpk {Heterogeneous nuclear 
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Figure 7 Venn diagram showing overlap between genes 
differentially expressed in all three conditions. Venn diagram 
showing annotated endogenous genes tliat are differentially 
expressed in all three conditions - current (VLMi), current + GFP 
(VLMg) and current + DmrtS (VLMd) - and the fold change range 
for each condition. 

V J 

ribonucleoprotein K), Chdl {Chromodomain-helicase- 
DNA-binding protein 1) and Txnll {Thioredoxin-like 
protein 1). Hnrnpk is an RNA binding protein known 
to associate with pre-mRNAs in the nucleus and influ- 
ences RNA processing, metaboUsm and transport It is 
involved in maintaining cellular ATP levels in stress 
conditions, possibly by protecting its target mRNAs 
[38]. Chdl is a chromatin remodelling enzyme that 
regulates transcription [39], and Txnll is a redox sen- 
sor [40]. 

Of the 21 genes aflected by current, 7 are also afl'ected 
by current + GFP and 12 are also affected by current + 
DmrtS (Figure 7). This indicates that the changes identi- 
fied in response to current are reproducible, but that 
the addition of DNA does alter the tissue's response to 
the current. Of the 111 genes affected by current + 
GFP, only 19 are affected by exogenous expression of 
DmrtS (Figure 7). This again indicates that the genes 
responding to the experimental conditions are affected 
in a combinatorial nature. 

Conclusions 

We have investigated the effects of both the current 
required for in ovo electroporation and the effect of 
exogenous DNA expression resulting from in ovo elec- 
troporation. The results of microarray analysis show 
that current and current + GFP have a small effect on 
endogenous gene expression, without affecting regional 
identity. This study demonstrates the validity of the in 
ovo electroporation technique to study gene function 
and expression within the developing embryo. 



Materials and methods 

In ovo electroporation 

HH stlO embryos were windowed and electrodes place 
either side of the developing head (CUY610P1.5-1, Soni- 
del Ltd (Dublin, Ireland). A solution containing 7 (ig/(il 
DNA construct pCAp-IRES-GFP or pCAp-IRES-DmrtS, 
2% polyvinyl alcohol, 0.05% Fast Green in water was 
injected into the developing midbrain. Five 12 V square 
wave pulses of 50 ms duration with an interval of 100 
ms were applied across the electrodes using an ECM830 
Electro-S Square Porator (BTX Inc. (Hawthorne, NY, 
USA). Following exposure to current, embryos were 
incubated at 37°C for 24 hours; embryos were collected 
at HH stl6. 

Tissue preparation 

Collected embryos were washed in PBS and transferred 
to clean PBS for dissection of the midbrain region, as 
marked by red Unes in Figure lA. Dissected tissue was 
transferred to 2 U/ml dispase (Sigma (St Louis, MO, 
USA) for 10 to 15 minutes. Mechanical dissection was 
used to remove the neural tube from the mesoderm. 
VLM was dissected from wild-type embryos (VLM), 
embryos that had experienced electric current (VLMi), 
embryos exposed to current and GFP (VLMg) and 
embryos exposed to current and the regulatory gene 
DmrtS (VLMd); Figure lA(ii-iv) shows schematics of the 
region isolated. Only embryos showing high expression 
of GFP were used; Figure lA(i) shows an in situ of an 
electroporated embryo, and Figure lA(ii-iv) show sec- 
tions of electroporated embryos, highlighting the elec- 
troporation efficiency and region that was dissected 
from the embryos. 

For each microarray, six VLM, six VLMi, six VLMg or 
six VLMd tissue samples were collected and pooled. 
Total RNA was isolated from each pool (Absolutely 
RNA Miniprep Kit, Stratagene (La JoUa, CA, USA). 
RNA was sent to UK Bioinformatics for processing and 
analysis. (London, UK) The procedure used was as fol- 
lows. cDNA was biotin labelled and samples were hybri- 
dized to Affymetrix GeneChip Chicken Genome Array 
according to the manufacturer's protocol. For each con- 
dition three sets of pooled samples were collected and 
hybridised to three arrays to obtain three biological 
replicates for each condition. 

Bioinformatics 

Probe levels were calculated from raw data using the 
MAS5 algorithm embedded into the GCOS suite (ver- 
sion 1.2; Affymetrix). Data were analysed using the 
GeneSpring package (version GX7; Agilent Technolo- 
gies, Wokingham, Berkshire, UK). The suitability of the 
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expression data sets for inclusion in the analysis and the 
overall relationship between biological replicates was 
assessed using quantile plots and PCA. Samples were 
first normalised to the 50th percentile (median) across 
the entire expression data set. Genes were then filtered 
to remove any genes absent in all arrays. To identify 
genes showing differential expression between the three 
samples, a one-way ANOVA with no multiple testing 
correction was carried out. A P-value of 0.05 and a two- 
fold change in gene expression was used to determine 
genes showing differential expression between the 
samples. 

Ingenuity network analysis 

The differentially expressed genes were analysed using 
IPA V8.0-2602 (Ingenuity Systems, Redwood City, CA, 
USA). The Ingenuity Knowledge Base is the largest 
knowledge base of its kind, with millions of findings 
curated from the full text literature. For analysis, the 
Affymetrix Chicken Genome Array gene list was used as 
a reference set. All data sources, all species, and all tis- 
sues and cell lines were used for the analysis. IPA uses a 
Fisher s exact test to determine which toxicity pathways 
and biological functions were significantly enriched 
within the set of genes showing differential expression 
compared to the entire list of genes represented on the 
array. 

Gene Ontology Tree Machine 

GO provides a controlled vocabulary of about 20,000 
terms in three independent hierarchies for cellular com- 
ponents, molecular functions and biological processes. 
The genes showing differential expression upon expo- 
sure to current or current + GFP were converted into 
mouse annotation and input into the GOTM program. 
As a reference set the entire mouse genome was chosen. 
The GOTM web-based tool carries out statistical analy- 
sis to identif)^ enriched GO categories for the input gene 
sets and generates a GO tree to visualize GO terms that 
are enriched in the input gene list. Hypergeometric test 
was used to select enriched GO terms for each cluster 
compared to the GO terms of the entire mouse genome. 
A GO category was considered as enriched if the P- 
value was <0.01. 

In situ hybridisation 

Whole mounts were exposed to digoxigenin-tagged anti- 
sense RNA mouse DmrtS probe overnight at 70°C for 
hybridisation with the exogenous DmrtS following the 
protocol described in [41]. 

Abbreviations 

CHO-Kl: Chinese hamster ovary epithelial cells; DmrtS: Double-sex and 
mab3-related transcription factor like 5; GFP: green fluorescent protein; GO: 



Gene Ontology; GOTM: Gene Ontology Tree Machine; HEK293: human 
embryonic kidney cells; HH: Hamburger-Hamilton; IPA: Ingenuity Pathway 
Analysis; IRES: internal ribosomal entry site; NIH3T3: mouse fibroblast cells; 
PBS: phosphate-buffered saline; PCA: principal component analysis; st: stage; 
VLM: ventral lateral midbrain; VLMd: ventral lateral midbrain exposed to 
current and Dmn5] VLMg: ventral lateral midbrain exposed to current and 
GFP] VLMi: ventral lateral midbrain exposed to current. 
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